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Previously observed non-Arrhenius behavior in fast ion conducting glasses [Phys. Rev. Lett. 76, 
70 (1996)] occurs at temperatures near the glass transition temperature, Tg, and is attributed to 
changes in the ion mobility due to ion trapping mechanisms that diminish the conductivity and 
result in a decreasing conductivity with increasing temperature. It is intuitive that disorder in 
glass will also result in a distribution of the activation energies (DAE) for ion conduction, which 
should increase the conductivity with increasing temperature, yet this has not been identified in the 
literature. In this paper, a series of high precision ionic conductivity measurements are reported for 
0.5Na2S + 0.5[a::GeS2 + (1 — a;) PS 5 ] glasses with compositions ranging from < a; < 1. The impact 
of the cation site disorder on the activation energy is identified and explained using a DAE model. 
The absence of the non-Arrhenius behavior in other glasses is explained and it is predicted which 
glasses are expected to accentuate the DAE effect on the ionic conductivity. 
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It is universally accepted that the ionic conductivity of 
a simple sohd electrolyte should obey an Arrhenius rela- 
tion. The temperature dependence of the conductivity is 
controlled by the activation energy, which is physically 
interpreted as the energy barrier that must be overcome 
for an ionic charge carrier to jump to an adjacent site. 
This simple relationship is expected in crystals where the 
activation energy is a singular value. In highly disor- 
dered systems, such as ionic glasses, the ion conduction 
processes are expected to arise from jumps over a distri- 
bution of such energy barriers, therefore, a distribution 
of activation energies (DAE) is anticipated. The atomic 
level origin of the DAE arises from the structural disorder 
in the glass, at both short and long length scales, which 
results in wide and continuous distributions of bond dis- 
tances, bond angles, and in some cases coordination num- 
bers. This DAE should result in the ionic conductivity 
having a non-Arrhenius temperature dependence, where 
the low energy barriers are crossed at low temperatures 
and the high energy barriers are crossed at high tem- 
peratures. Astonishingly, however, the vast majority of 
published studies of ion-conducting glasses report an Ar- 
rhenius relationship with a single activation energy. Until 
now it remained unknown why the chemical and struc- 
tural complexity that is intrinsic to the glassy state does 
not result in a dramatic deviation from simple Arrhe- 
nius behavior. In this letter we resolve this long standing 
question. 

There are a limited number of studies observing non- 
Arrhenius conductivity in glasses. One class of non- 
Arrhenius behavior is exemplified by the measurements 
of Kincs and Martin on silver halide doped sulfide glasses 
[H. At low temperatures, the conductivity maintains a 
linear Arrhenius behavior, but at higher temperatures 
it develops a distinctly negative curvature. The models 



developed to explain this behavior have all focused on 
site-hopping diffusion and mobility [INBl- Malki et al. 
identify a rigid-to-floppy transition, which depends on 
both the composition and the temperature to delineate 
the high-temperature, non-Arrhenius behavior from the 
low temperature regime In essence, this form of 

non-Arrhenius behavior is caused by a temperature de- 
pendent mobility due to an ion trapping mechanism that 
reduces the mobility at high temperatures 0-0] • 

A second class of non-Arrhenius behavior is character- 
ized by a positive curvature of the conductivity across 
the entire temperature region. Archetypal examples of 
this behavior are found in Namikawa's measurements of 
mixed-alkali glasses and Murugavel et al. for mul- 
ticomponent phosphosilicate glasses Namikawa at- 



tributes this behavior to the existence of multiple charge 
carrying species and uses two straight lines to approxi- 
mate the curved behavior, with each slope indicating the 
activation energy of one of the species. While it is pos- 
sible to attribute the positive-curvature non-Arrhenius 
behavior to the activation energies of two charge car- 
riers; the same non-Arrhenius behavior is observed in 
glasses with single charge carriers: Murugavel and Rol- 
ing for a;Na20 -I- (1 — a;)B203 and Imre et al. for 
0.2Na2O + O.8B2O3 and 0.2Rb2O + O.8B2O3 [H. 

In this letter we report on a series of high preci- 
sion measurements on mobile Na ion containing sul- 
fide glasses. The data demonstrates the same positive- 
curvature non-Arrhenius behavior described in refer- 
ences 13- 13| • Following the lead of our earlier work [B-Q , 
we develop a DAE model and demonstrate an excellent 
agreement between our model and the data. The results 
of this model allow us to explain why such non-Arrhenius 
behavior is not commonly seen in other glasses and how 
to modify a glass to enhance the influence of the DAE. 
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Finally, DAEs have been used to describe ultrasonic at- 
tenuation (UA) measurements of glass where the UA loss 
peaks were observed to be wider than could be described 
by a single activation energy (l4. 15|. 

Mixtures with compositions 0.5Na2S + 0.5[a;GeS2 + 
(1 — x) PSs], where X ranges from 0.0 to 1.0, were pre- 
pared from' O.SNaaS-h 0.5 GeSz andO.5Na2S-HO.5PS5. Ah 
syntheses were carried out in a high purity N2 glove box. 
The mixtures were melted in covered vitreous carbon cru- 
cibles inside a tube furnace at 730 °C for 3 minutes. The 
mixture was removed from the furnace and allowed to 
cool inside the crucible. The mass loss due to sublima- 
tion of reactants, typically less than one percent, was 
recorded to ensure compositional accuracy. Disk-shaped 
specimens were then prepared by remelting the glass and 
pouring it into disk-shaped brass molds that were held 30 
°C below the glass transition temperature. The glasses 
were transparent and had yellow to amber colors. After 
annealing for 30 minutes, the disks were cooled to room 
temperature at 1 °C/min. The disks were polished and 
circular gold electrodes were sputtered onto both faces 
of the glass. A mask was used to assure that the elec- 
trodes were registered on both sides of the disk. Using 
a Novocontrol Technologies Concept 80 impedance spec- 
trometer, the complex impedance spectra were measured 
from 0.1 Hz to 3 MHz and temperatures from -50 °C to 
150 °C, where -50°C is the lower limit of the components 
of the custom cell used for O2 and H2O sensitive materi- 
als and 150°C is below the annealing temperature for all 
samples in this study. The temperature was held within 
±0.5 °C of the nominal set point for three minutes prior 
to data collection to stabilize the temperature. The di- 
rect current (DC) conductivity was found by fitting the 
complex impedance arc and using the disk's thickness 
and electrode area. 

The ionic conductivity of glass can be expressed as the 
Nernst-Einstein relation, a modified functional form of 
which is, 



from 
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where R is the gas constant and the prefactor, (Tq, con- 
tains the ion charge, number density, and diffusivity. The 
activation energy, AEa, is an energy barrier that the 
charge carrying ions must overcome to conduct through 
the glass network. In conductivity studies, it is common 
practice to use the measurements to determine the values 
of (Jo and AEa as constant material properties. However, 
this assumes that both parameters are singular constant 
values, whereas the disordered short and intermediate 
range chemical structures seen by mobile ions in glasses 
should result in AEa being a DAE associated with these 
local variations. The simplest correction to this assump- 
tion is to replace AEa in equation [T] with a temperature 
dependent expectation energy, (A_B). that can be found 
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where P {AE, T) is the probability distribution function 
of the activation energies of the mobile ions at temper- 
ature T. This probability distribution introduces a tem- 
perature dependence to the expected energy and is writ- 
ten using the Boltzmann relation, 
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The function g (AE), in equation |3l is the temperature- 
independent probability distribution function of the ac- 
tivation energies, DAE, of all the ions in the glass. In 
this method, it is assumed that the energy landscape 
seen by the mobile cations is fixed upon quenching and 
with increasing temperature the number of these fixed en- 
ergy barriers being overcome increases according to the 
Boltzmann distribution. For a given g (AE) , equation [3] 
is integrated to determine the probability distribution 
function of the mobile ions as a function of tempera- 
ture, P{AE,T). Subsequently, equation [2] is used to de- 
termine the temperature dependent expectation energy, 
( AE') , which then is used in place of AEa in equation [1] 
The measured ionic conductivity of a typical glass in 
this series is plotted in FIG. [Ija). On first inspection, 
the data appears to be Arrhenius; however, closer evalu- 
ation reveals a positive, upward curvature. This non- 
Arrhenius nature is more evident when the data are 
subtracted from a straight line drawn between the low- 
est and highest temperature data points, as shown in 
FIG. [TJb). The deviation from linearity is distinctly 
quadratic and this has been observed by Imre et al. for 
0.2[a;Na2O + (1 - a-)Rb20] + O.8B2O3 glasses which is 
plotted in FIG. 8 of reference [l^l. The non- Arrhenius be- 
havior demonstrated in FIG. [T] is systematically present 
for all 10 glass specimens with compositions ranging from 
a; = to 1. 

For these glasses, the DAE barriers, g{AE), can be 
taken to be Gaussian, with a mean value of AEq and 
standard deviation of S. 0-Q This results in equation [5] 
being expressed. 
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which can be simplified to 
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In this way, the temperature-dependent (AE) can be 
written in terms of the Gaussian parameters, AEq and 
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FIG, 1, The measured ionic conductivity for the glass com- 
position 0,5Na2S + 0,5[0.7GeS2 + O.SPSs], Frame (a) shows 
the conductivity plotted with a best-fit Arrhenius line, where 
the experimental error is smaller than the symbols. The data 
in frame (b) is the difference between the conductivity data, 
from (a) minus a point from a line arbitrarily drawn between 
the highest and lowest data points. The solid line in (b) is 
the match between the Arrhenius fit in (a) and the arbitrary 
straight line. The dashed line is the result of the DAE model. 



6. To write the simplified expression in equation |3] it is 
required that P {AE, T) is zero at T = 0. This approx- 
imation breaks down in the hmit of smah T and large 
6/Eq. The glasses and temperatures used in this study, 
and in most others, are well within the operating limit of 
this approximation. The deviation between the full and 
simplified expressions is shown as dotted lines in FIG. 
HJb). This elearly demonstrates that for T > 50 K equa- 
tionlholds for S/AEq = 0.067. 

Using this temperature dependent expectation energy 
in place of the constant activation energy in equation [l] 
the natural log of aT becomes quadratic in 1/T. For 
each composition of glass, the parameters AEq, S, and 
ctq were determined by replacing AEa in equation [T] with 
(AE) from equation 3] and best-fitting this expression to 
the experimental data. Results of this fitting are shown 



in FIG. HJa). See Supplemental Material at [URL will 
be inserted by publisher] for a table of DAE parameter 
values for all of the glasses. As can be seen by the fitted 
curve in FIG.HJb) the model matches the data extremely 
well. 

To understand the impact of the DAE on the thermal 
behavior of the ionic conductivity, {AE) is parametri- 
cally plotted for different S/AEq ratios as a function of 
temperature in FIG. IH^b). The experimental tempera- 
ture range of this study lies between the dashed lines. 
For this family of curves, it is observed that the temper- 
ature dependence is greater at lower temperatures and 
larger S/AEq ratios; hence, the non- Arrhenius behavior 
will be more pronounced in material systems with DAEs 
that have a small mean activation energy, AEq, and a 
large width, 5. Glass compositions with large AEq val- 
ues such as the most commonly studied silicate, borate, 
and phosphate oxide glasses are not expected to exhibit 
the non- Arrhenius behavior shown in FIG. [U even if the 
6 is equivalent to the values found in analogous non-oxide 
glass systems, such as the glasses in this study. The non- 
Arrhenius behavior will be enhanced in highly conductive 
materials with low AEq values. For example, our model 
predicts that glasses that have Li as the charge carrying 
species will demonstrate a stronger DAE non- Arrhenius 
behavior than the other alkali ions because the small size 
of the Li ions results in a lower mean activation energy, 
AEq. An example of this can be found in FIG. 4 of refer- 
ence [3l where the ionic conductivity of 36Li2S ■ 18P -468 
is shown to have a highly non- Arrhenius conductivity be- 
low the Tg. 

At lower temperatures, there is a stronger temperature 
dependence of the expectation energy because only the 
low energy tail of g {AE) contributes to P{AE,T) via 
equation [31 As the temperature increases, the expecta- 
tion energy asymptotes to AEq because more of the en- 
ergies in the g {AE) distribution begin to participate; at 
high temperatures the expectation energy becomes con- 
stant and the measured conductivities will appear Ar- 
rhenius. This implies that low-temperature experiments 
will accentuate the observation of the DAE. The exact 
experimental range needed is dependent on the material. 

The non- Arrhenius behavior also may be enhanced 
by the processing methods used to prepare and the 
chemistries of the glass, both of which can be used to 
increase S. Rapid quenching of a glass melt will freeze 
in structures with a higher Active temperature and, thus, 
a higher configurational entropy. RF-sputtered glassy 
thin films may also lead to a wider distribution of atomic 
level disorder. Additional disorder, and possibly coordi- 
native defects, can be introduced by mechanical means 
such as high energy ion implantation or planetary milling. 
Glasses with greater chemical complexity will also lead 
to increased 6. A distribution of anionic sites contributes 
to the DAE, and quench rates are known to strongly in- 
fluence speciation in glasses |17| . Further, the use of mul- 
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FIG. 2. Frame (a) shows the room temperature ionic con- 
ductivity and DAE model parameters, AEo and S, for glasses 
with composition 0.5Na2S + 0.5[a;GeS2 + (1 — a;) PS 5] as a 
function of composition, where the lines are guides for the 
eyes and error bars are smaller than the symbols. Frame (b) 
shows the expectation energy as a function of temperature for 
a variety of 5/ AEo choices. This family of curves asymptotes 
toward AEq with increasing temperature. The rate of con- 
vergence depends on the value of 5/ AEq. The dashed lines 
bracket the temperature range for the glass specimens studied 
here. The solid lines show the exact integral, from equation[51 
and the dotted lines show the expectation energy from equa- 
tion [l] The data plotted in FIG. [T] corresponds to a; = 0.7 
which has 6/ AEq = 0.067. 



tiple charge carriers will increase the width of the DAE, 
and the quadratic deviation from Arrhenius behavior ob- 
served in reference 
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indicates that the existence of 
multiple charge carriers is not a simple matter of having 
two sin gula r valued activation energies as Namikawa pos- 
tulates [10[. If only two activation energies exist, then the 
deviation from linearity should have a simple "V" shape 
instead of being quadratic. 

Glasses with more than one glass former, such as the 
mixed glass former glasses in this study, should also have 
a wider DAE. FIG. HJa) shows the compositional de- 
pendence of the ionic conductivity, AEq, and 6 in the 
0.5Na2S -I- 0.5[a;GeS2 + {1 - x) PSs] glasses. As is ex- 
pected, the composition with the largest AEq, x = 0.7, 



has the lowest conductivity. The structural origins of the 
compositional dependence of AEq have so far not been 
fully identified, but the compositional dependence of 5 
can be linked to the distribution of local structures. See 
Supplemental Material at [URL will be inserted by pub- 
lisher] for the relative abundance of short-range struc- 
tures in the 0.5Na2S + 0.5[a;GeS2 + {1 - x) PSs ] glasses. 
Structures are labeled Q^, where Q denotes the glass- 
forming cation (P or Ge) and z denotes the number 
of bridging sulfurs associated with the structural unit. 
When a; = 0, the dominating structural unit is the P""^ 
group, nominally Na2PS 7 . Addition of Gc^ groups, nom- 
inally Na2GeS3, leads to formation of P'^ and Ge'^ groups 
arising from a disproportionation reaction 



+ Ge^ P" + Ge'. 



(5) 



The consequence of this proposed reaction is that Na 
ions experience a wider distribution of chemical envi- 
ronments. Qualitatively, ternary compositions appear 
to have a wider distribution of structures where no 
group is dominant. For example, the x = 0.3 composition 
is comprised of 20.6% P\ 16.4% P^^, 32.1% po, 29.6% 
Ge'^, and has the largest 5 value of all of the glasses. 

In this letter, we have examined the impact that the 
DAE has on the ionic conductivities in glasses by per- 
forming a series of high precision measurements on a Na 
containing mixed-glass former sulfide glasses and con- 
structing a model that shows how the DAE will affect 
the Nernst-Einstein relation. Using this model, we have 
addressed the question of why the positive-curvature non- 
Arrhenius conductivity associated with the DAE has only 
been identified very infrequently. It is significant that of 
the many thousands of ionic glass systems studied only 
a small number have reported this behavior. The most 
studied glasses are the more poorly conducting, yet most 
easily prepared, oxide glasses whose large AEq and small 
5 values produce essentially an Arrhenius conductivity. 
It is only in materials such as the more highly conduct- 
ing non-oxide glasses, especially the complex mixed-glass 
former glasses, that the AEq is sufiiciently low and the 
5 is large enough that the non-Arrhenius conductivity 
is apparent. This study is the first to understand and 
report the underlying cause of this apparent unknown 
behavior. 
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